A new class of organic molecular magnets 
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Using detailed first principles calculations, we have found a new class of stable organic molecular 
magnets based on zwitterionic molecules possessing donor, vr bridge, and acceptor groups. These 
molecules are organic molecules containing only C, H and N. The quantum mechanical nature of 
the magnetic properties originates from the conjugated tt bridge (involving only p electrons) where 
the exchange interactions between electron spin are relatively strong and local and are independent 
of the length of the tt bridge, enabling the easy construction of a molecular magnetic device with 
specified length. 



PACS numbers: 75.50.Xx 

Metal-free molecular magnets offer intriguing benefits 
in terms of cost, weight, and diversity in magnetic proper- 
ties and ease in processing and synthesis. However most 
traditional molecular-based magnets are crystals com- 
posed of molecular units that contain a variety of metals 
or radicals [T] that have either a low Curie temperature 
(Tc) or are unstable due to the processing techniques 
used to generate them. In the present paper we show that 
a particular class of organic molecules, containing no rad- 
icals or metals, can have a stable magnetic ground-state 
and thus, should have a high Neel temperature (Tjy)- 
Moreover, by simply substituting a N with a C atom 
on the 6-membered ring comprising the betaine molecule 
(see figure [ij the magnetic state can be changed from an- 
tiferrimagntic (Afm) to ferrimagnetic (fm) and it should 
have a high Tp. Magnetic materials that are purely or- 
ganic molecules are uncommon since the p-electrons in 
carbon based materials generally participate in covalent 
bonds [2] . Our findings indicate that betaine derivatives 
promise to offer a new class of molecular magnetic materi- 
als with diverse potential applications including organic 
spintronics, information storage and nanoscale sensors. 

A limited number of molecular-based magnets have 
been reported [T] and these are typically composed of or- 
ganic radicals [3] or mixed coordination compounds con- 
taining bridging organic radicals, Prussian Blue type 
compounds, or charge transfer complexes [3|. For the 
relatively few organic materials that possess magnetic 
properties, the Curie temperature (Tc) is generally very 
low (<36 K) [5][S][B] while those with high Tc typi- 
cally contain some metallic atomsfT^ f8^ [9] |2] with the or- 
ganic molecular part of the crystal acting as exchange 
pathways. Other allotropes of carbon such as carbon 
nanofoams[TU| and boron-nitride nanotubesjTI] have also 
been reported to exhibit magnetic properties, however 
either the magnetic properties are very short-lived or 
they depend on doping and defects to achieve a mag- 
netic state where the stability and strength of the mag- 




FIG. 1: A zwitterionic betaine molcule with an imidazol 
donor (D) and a pyridine acceptor (A) separated by alternat- 
ing C=C bridge. 



netism is low since the achieveable doping/defect con- 
centration is quite low[12 . Organic free radicals are 
natural candidates for magnetic materials, but very 
few are stable enough to be used as magnetic devices. 

While it is clear that some progress towards an or- 
ganic magnetic material has been made (in particular 
for graphitic systems), room temperature, air stable and 
synthetically viable molecules are still lacking. Here we 
report on a zwitterionic molecular system that is com- 
posed entirely of C, N and H and has a magnetic ground- 
state. The betaine derivatives are zwitterions with a tt- 
conjugated structure whose variable tt bridge length sep- 
arates the donor and acceptor groups [TS]. 7r-conjugated 
organic molecules are in general good candidates for ex- 
hibiting magnetic order, [H] although as discussed above, 
only a relatively few successful examples have been re- 
ported and those have very low Tc. The betaine deriva- 
tives are unique because they also have an intrinsic donor 
and TT acceptor at different ends of the molecule. As 
will be shown in this paper, the magnetic properties 
of these systems are independent of the bridge length. 

The molecular structure for the betaine derivatives is 
shown in Figure [l] . These systems are zwitterions, con- 
taining a pyridine acceptor (cation) and imidazole donor 
(anion) separated by a variable length 7r-conjugated 
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(C2H2)„ bridge. In this paper, we utilize a combina- 
tion of ab initio techniques for performing fully relaxed 
total energy calculations using Hartree-Fock (HF), com- 
plete active space self consistent field (CASSCF), and 
density functional theory (DFT) within the framework 
of the NWChem package [TS], to determine the intrinsic 
character of the ground state for several betaine deriva- 
tives. The standard Pople basis set, 6-311G** was used 
for all calculations, yielding similar results. A hybrid ap- 
proach was chosen since it is known that using standard 
density functional theory methods for purely organic 
TT conjugated systems often fail due to the delocaliza- 
tion of the system [16], and 2) an accurate description of 
electron exchange is necessary for determining magnetic 
properties. We therefore only report the DFT results 
obtained using 100% HFX with Perdew-Zunger (PZ) 
correlation 1 17j. We note that one can use either PZ or 
Perdew-Burke-Ernzerhof (PBE)[T8] correlation without 
affecting the magnetic properties since the mechanism 
responsible for antiferrimagnetism or ferrimagnetism is 
purely local. To determine the effects of non-dynamical 
correlation on the magnetic properties, we used CASSCF 
and the 6-31 IG** basis set with an active space of up to 
14 active electrons and 14 orbitals. The CASSCF re- 
sults are in good agreement with those obtained from 
HFX-DFT indicating that the HFX-DFT approach is 
adequate in describing the correlation in this system. 

For the ferromagnetic (FM) state several multiplicities 
were considered and full geometry optimization was per- 
formed where the triplet state was found to be the lowest 
in energy for the chain lengths considered. We note that 
the initial antiferromagnetic state relaxes to an antiferri- 
magnetic state (Afm) , and this is the lowest energy state 
for all bridge lengths examined. Similarly, the initial FM 
state relaxes to a ferrimagnetic state (fm) for all bridge 
lengths considered. The results of our investigation 
demonstrate that the ground state is magnetic, imply- 
ing that this class of molecules are excellent candidates 
for the development of new organic magnetic systems. 

Restricted open-shell Hartree-Fock (ROHF) calcula- 
tions for the fm type states were used to rule out potential 
spin-contamination problems. Higher multiplicity calcu- 
lations for the fm state were also performed to determine 
the lowest energy spin state. For the fm type calculations 
ROHF, UHF and HFX all obtained the triplet state as 
the ground state spin configuration and no significant 
effects associated with spin contamination were found. 

Table 1 lists the results obtained from geometry op- 
timized calculations using UHF and HFX-PZ. Both 
yield the same spin configuration ordering with the 
lowest energy state being antiferrimagnetic (Afm), fol- 
lowed by a ferrimagnetic state (fm) with the highest 
energy state being non-magnetic (NM). The Afm-fm 
energy differences {AF,Afm-fm) listed in table 1 indi- 
cate that the Afm ground-state spin structure is ex- 
ceptionally stable, exhibiting an equivalent temperature 



TABLE I: Total energy difference for fully relaxed UHF and 
HFX-PZ (labeled PZ) simulations for NM, fm and Afm states 
along with the Afm-NM and Afm-NM and fm-NM total enery 
differences for UHF and PZ respectively. Atomic units are 
used and n in Betn, defines the bridge length. 





UHF 


AEAfm-fm 


PZ 


AEAfm-fm AEfm-NM 


n = 2 










nm 


-624.064 




-630.792 




fm 


-624.105 




-630.811 


-0.0187 


Afm 


-624.117 


-0.0118 


-630.824 


-0.0127 


n = 4 










nm 


-777.843 




-786.352 




fm 


-777.919 




-786.387 


-0.0350 


Afm 


-777.921 


-0.0114 


-786.399 


-0.0127 


n = 5 










nm 


-854,734 




-864.1333 




fm 


-854.812 




-864.1744 


-0.0411 


Afm 


-854.823 


-0.0113 


-864.187 


-0.01234 


n = 6 










nm 


-931 626 




-941 916 




fm 


-931.713 




-941.963 


-0.0470 


Afm 


-931.726 


-0.0130 


-941.975 


-0.0124 


n = 8 










nm 


-1085.410 




-1097.481 




fm 


-1085.519 




-1097.539 


-0.0583 


Af 


-1085.530 


-0.0113 


-1097.551 


-0.0124 


n = 10 










nm 


-1239.195 




-1253.046 




fm 


-1239.322 




-1253.116 


-0.0693 


Afm 


-1239.335 


-0.0130 


-1253.128 


-0.0124 



that is well above room temperature. Furthermore, the 
A^Afm-fm is relatively constant as compared to ei- 
ther the AFijifjn-NM or AEf^-NM energy differences. 
This is about an order of magnitude larger than that 
typically reported for graphene nanoribbons [T^ |20j . 

Figure [2] shows the results for 3 different betaine 
molecules in the Afm and fm states with bridge lengths 
of 4, 5, and 10 (see figure [l] for a description of a 
chain length). The magnetic moments for each C and 
N site and the total magnetic moments for the chain, 
5 and 6 member rings are also displayed (the magnetic 
moments were calculated using a MuUiken population 
analysis (21]). For both magnetic states, the magnetic 
moments of the rings and chains as a function of in- 
creasing chain lengths are constant. Interestingly, in- 
dependent of the magnetic state the chains and rings 
exhibit an Afm and fm type ordering respectively. For 
the Afm state, the contributions from the rings approxi- 
mately cancel each other out, while for the fm state, only 
the magnetic moment of the 6-membered ring changes 
and it varies by ~ 1/iB. These results are the same 
for all chain lengths investigated in this work. Thus, 
it appears that the magnetism of the system can be 
viewed as consisting of three nearly independent com- 
ponents, two ring structures and an unsaturated chain. 

The above results imply that the chain length compo- 
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4-Afnn 5-Afm 10-Afm FIG. 4: AEA/m-iVM vs C2H2 vr bridge length. 



FIG. 2: The magnetic moments for C and N along with 
the 5,6 membered ring and C2H2 it bridge magnetic moment 
contributions for chain lengths of 4,5 and 10. 




FIG. 3: Figure 3a shows the delocalization of the charge den- 
sity. Figure 3b displays the spin density showing the strong 
Heisenberg character. 

nent acts as a strongly correlated Afm insulator. This can 
be seen in the spin density plot shown in figure [3]3. Note, 
that the spin density along the chain contains significant 
localized Heisenberg spin type character and does not 
contain any spin dimer like qualities, [H] [23] [21] while the 
charge density in figure [3]i exhibits strong 7r-conjugation 
throughout the molecule. The only fm like coupling oc- 
curs in the 5 member ring, which by symmetry must 
occur. To investigate the magnetic interaction further, a 
plot of the Afm-NM energy difference (AE^y,„_jvM) vs 
C2H2 chain length is shown in figure |4] For chain lengths 
of 4 to 10, AEAfm-NM decreases linearly with increasing 
chain length, indicating a strong Heisenberg like behav- 
ior with the change in energy being oc —JJ2{ij) ' ^j^ 
where the summation is restricted to nearest-neighbors. 
From the spin density and the magnetic moments in fig- 
ure [2] it is clear that the two rings only interact locally 
with the first couple of carbon atoms in the chains, which 
is why there is a slight change in linearity for the chain 
length of 2. The linear behavior of AEiAfm-NM with in- 
creasing chain length also indicates it is stable and that 
the chain acts to insulate the two ring structures. This 
type of screening would account for the change of 1 /iB 



on the 6-membered ring in the fm state. This can be 
seen from figure [2] where the moments of the two rings 
are essentially constant with increasing chain length. We 
also decomposed the eigenvalues (orbital and spin) onto 
atomic sites for the various chain lengths. The energy for 
the highest occupied states for the different chain lengths 
is -0.07 a.u. and essentially constant with increasing 
chain length. Furthermore, these states are of Pz char- 
acter and a strong p^ interaction occurs between neigh- 
boring spin-up and spin-down states within the chain. 
This feature taken along with figure 2 accurately illus- 
trates the strong localized Heisenberg spin character. It 
is also important to note that geometric defects in the 
chains, such as cis conformations, do not disrupt the Afm 
ordering (results are not shown). While substitutional 
defects that disrupt the tt conjugation on the bridge, 
such as saturation, will destroy the magnetic character. 

The Afm state in this system is different from Ander- 
sons theory of superexchange [25[ where two next-nearest 
neighbor transition metal atoms interact through a dia- 
magnetic anion. In addition, C and N atoms have their 
respective Px and Py orbitals shifted down in energy with 
respect to the Pz orbitals. On the other hand, the Pz 
spin orbitals are symmetry split, such that, the spin al- 
ternates between carbon sites on the chain. This order- 
ing combined with the local nature of the spins lowers 
the Coulomb energy between sites. The 7r-conjugation in 
the rings is strictly local to the individual rings and is 
maintained through primarily nearest-neighbor interac- 
tions with the chain. Interactions with neighbors further 
down the chain decrease almost exponentially, produc- 
ing a short-range effective screening, which again high- 
lights the role of the chains to insulate any interac- 
tions between the rings, even for short chain lengths. 

Additional support for the existence of a magnetic 
ground state was obtained from CASSCF calculations 
using various sized active spaces, ranging from 10 elec- 
trons and 10 orbitals to 14 electrons and 14 orbitals. 
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From these calculations it is clear that the natural or- 
bitals show significant partial occupancies, indicating a 
noteworthy multi-configurational character of the singlet 
ground state. For Bet6 with an active space of (12,12) the 
unpaired occupancies are: 1.9448, 1.9332, 1.9242, 1.8799, 
1.8570, 1.7075, 0.2933, 0.1525, 0.1202, 0.0739, 0.0523, 
0.0613. The unpaired electron density is localized on 
each C site across the C=C bridge and on the donor and 
acceptor rings of the molecule. For the larger betaines 
7-10, a similar result is found. Larger active space cal- 
culations did not change the fact that there is an open 
shell singlet ground state. The triplet state is higher in 
energy than the Afm singlet but below the non-magnetic 
singlet state. These results are consistent with those 
given in table 1 for the HFX-DFT and UHF calculations. 

In this paper we have presented a detailed examination 
of the fundamental electronic structure of a new class 
of zwitterionic molecules based on betaine derivatives. 
The results show that these compounds have a surpris- 
ingly strong local Heisenberg type magnetic character 
that originates from the conjugated bridge in coopera- 
tion with the donor-acceptor character of the molecule. 
The localized behavior found in this system is compat- 
ible with the TT-conjugation of the components making 
up the system, thus maintaining the resonant behavior 
in the systems charge density aiding in its overall sta- 
bilization. The Afm state further stabilizes the system 
because the Coulomb energy is diminished by alternat- 
ing the spins throughout most of the system. Thus, the 
stabilization of the system is achieved by the resonating 
nature of the charge density (kinetic energy) and the lo- 
calization of the Afm state that results in a lowering of 
the Coulomb energy. The Afm type ordering found on 
the conjugated bridge is critical in establishing the short- 
range local interactions with the donor and acceptor ring 
structures and the resulting lack of long-range interaction 
between them. These features are responsible for the sta- 
ble magnetic configuration as a function of bridge length 
and could possibly enable the development of a variety 
of magnetic devices. Furthermore the magnetic charac- 
ter is not disrupted by conformational disorder. Stable 
room temperature molecular magnetic materials that are 
purely organic are rather rare and the betaine derivatives 
promise to offer a new class of molecular magnetic materi- 
als with diverse potential applications including organic 
spintronics, information storage and nanoscale sensors. 
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